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ABSTRACT 

An attractive scenario for producing Type la supernovae (SNe la) is a double detonation, where detonation 
of an accreted helium layer triggers ignition of a C/O core. Whether or not such a mechanism can explain some 
or most SNe la depends on the properties of the helium burning, which in turn is set by the composition of the 
surface material. Using a combination of semi-analytic and simple numerical models, I explore when turbulent 
mixing due to hydrodynamic instabilities during the accretion process can mix C/O core material up into the 
accreted helium. Mixing is strongest at high accretion rates, large white dwarf (WD) masses, and slow spin 
rates. The mixing would result in subsequent helium burning that better matches the observed properties of 
SNe la. In some cases, there is considerable mixing that can lead to more than 50% C/O in the accreted layer 
at the time of ignition. These results will hopefully motivate future theoretical studies of such strongly mixed 
conditions. Mixing also has implications for other types of WD surface explosions, including the so-called .la 
supernovae, the calcium-rich transients (if they arise from accreting WDs), and metal-enriched classical novae. 
Subject headings: accretion, accretion disks — hydrodynamics — white dwarfs 


1. INTRODUCTION 

The use of Type la supernovae (SNe la) as cosmic distance 
indicators (e.g., Riess et al. 1998; Perlmutter et al. 1999) has 
put increased emphasis on the theoretical uncertainties that 
still remain about these events. It is generally accepted 
that SNe la result from the unstable thermonuclear igni¬ 
tion of a degenerate C/O white dwarf (WD), but the spe¬ 
cific progenitor systems have not yet been identified. The 
two main candidates can roughly be divided into (1) sta¬ 
ble accretion from a non-degenerate binary companion un¬ 
til the Chandrasekhar limit is reached (single degenerates, 
Whelan & Iben 1973), (2) the merger of two C/O WDs (dou¬ 
ble degenerates, Iben & Tutukov 1984; Webbink 1984). An 
important problem is therefore to understand which scenar¬ 
ios explain the SNe la that we observe and whether any one 
channel is dominant. 

How these systems ignite is also an important outstand¬ 
ing question. One ignition scenario that could occur in both 
single and double degenerate systems is the “double detona¬ 
tion.” This happens when accretion and detonation of a he¬ 
lium shell on a C/O WD leads to a prompt detonation of the 
core (Woosley & Weaver 1994; Livne & Arnett 1995). Dou¬ 
ble detonations have a number of attractive features which has 
led to increased interest in this scenario in recent years. Ob- 
servationally, they are consistent with limits on shock cooling 
of the exploding WD (Piro et al. 2010; Bloom et al. 2012), 
the non-detection of a companion in pre-explosion imaging 
of nearby SNe la (Li et al. 2011), the lack of radio emission 
(Hancock et al. 2011; Horesh et al. 2012), the lack of hydro¬ 
gen emission in nebular spectra (Leonard 2007; Shappee et al. 
2013), a lack of a signature of ejecta interaction with a 
companion (Kasen 2010; Hayden et al. 2010), the missing 
companions in SNe la remnants (Schaefer & Pagnotta 2012), 
and the sub-Chandrasekhar ejecta masses inferred in many 
cases (Scalzo et al. 2014). On the theoretical side, the sub- 
Chandrasekhar masses of the C/O cores when they ignite 
can roughly explain the size of the width-luminosity rela¬ 
tion (Phillips 1993) observed for typical SNe la (Sim et al. 
2010) and population synthesis calculations argue that dou¬ 


ble detonations could explain the majority of normal SNe la 
(Ruiteretal. 2011, 2014). Substantial theoretical work has 
also focused on the details of the helium and/or C/O burning 
(Fink et al. 2007, 2010; Moll & Woosley 2013). 

The main issues with the double detonation scenario have 
been whether the helium ignition is robust (Holcomb et al. 
2013) and whether the colors and spectra due to the ashes left 
over from the surface helium burning are strongly inconsis¬ 
tent with normal SNe la (Woosley & Kasen 2011; Sim et al. 
2012). However, it was recently argued that both of these dif¬ 
ficulties may be overcome by if there is additional C/O mate¬ 
rial in the surface helium layer prior to when the detonation 
begins (Kromer et al. 2010; Shen & Moore 2014, the latter of 
which also explored a more detailed treatment of the nuclear 
network). The main argument is that the additional burning 
decreases the hotspot size and mass needed for triggering a 
detonation, which in turn leads to less iron-peak elements in 
the helium layer that adversely affect the spectra. In addition, 
the silicon and calcium left over from the surface burning may 
lead to high-velocity spectral features as observed in most 
SNe la (Mazzali et al. 2005; Tanaka et al. 2008; Blondin et al. 
2012; Childress et al. 2014; Maguire et al. 2014). This C/O 
can get created in the surface layer by the convective he¬ 
lium burning in the time prior to initiation of the detona¬ 
tion (Shen & Bildsten 2009) or mixed in from the C/O core 
by convective mixing associated with this same stage. Here 
I consider in more detail whether C/O material can also be 
mixed into the helium layer via turbulent mixing during the 
accretion process. 

Additional motivations for this work are the scenarios 
where the C/O core is not successfully ignited by the sur¬ 
face helium detonation (Shen & Bildsten 2014). In such 
cases, the low mass of the helium would lead to a faint and 
rapidly evolving transient (Bildsten et al. 2007; Shenetal. 
2010; Waldman et al. 2011), which may be related to var¬ 
ious recently observed rapid transients (Perets et al. 2010; 
Kasliwal et al. 2010, 2012; Foley et al. 2013; Inserraetal. 
2015). Whether or not a helium shell-burning scenario can 
explain a given observation will depend on the composition 
of the surface layers, which again can be affected by turbulent 
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mixing. 

The impact of spin on helium-accreting WDs has been the 
subject of previous theoretical work (Yoon & Langer 2004a; 
Yoon et al. 2004). These studies were mainly focused on the 
shear instabilities for providing the angular momentum trans¬ 
port, viscous heating, and material mixing. Since I show 
below that shear instability is likely subdominant to other 
processes, it is worth revisiting these models. Furthermore, 
new questions about the double detonation scenario and new 
classes of optical transients (both described above) mean that 
it is useful to address mixing in the context of these new po¬ 
tential applications. 

In Section 2, I summarize the picture of turbulent mixing 
and angular momentum transport that I utilize. In Section 
3, I discuss hydrodynamic mixing mechanisms and analytic 
scalings which demonstrate how it depends on accretion rate 
and spin. In Section 4, I compare with numerical models to 
demonstrate these main features. I conclude in Section 5 with 
a summary of my results and a discussion of future work. 

2. TRANSPORT AND MIXING FRAMEWORK 

Before moving into the main calculations, it is helpful to 
first explain the general strategy utilized by this work. Both 
the treatment of the angular momentum transport and the tur¬ 
bulent mixing are somewhat different than what is typically 
used, so it is helpful to explain why this approach is useful. 
Similar methods were used in Piro & Bildsten (2007), and 
thus readers interested in additional details should consult this 
work. 


2.1. Angular Momentum Transport 


the cases considered here and I instead focus on v as a turbu¬ 
lent diffusivity. 

Equation (3) can be significantly simplified by assuming 
steady-state transport of angular momentum (also see Piro 
2008). This is reasonable as long as the angular momentum 
diffusion timescale at a given depth 

ivisc = H 2 /v , (4) 

is less than the timescale to accrete down to that same depth 

f acc = 4 t rf? 2 X/M, (5) 

which is easily satisfied in the WD surface layers. The total 
time derivative then simplifies to d/dt « I4d v d/dz, where 
V^dv = — M/4nR 2 p is the advecting velocity of the fluid in 
an Eulerian frame. Making this substitution, integrating with 
respect to z, and taking the limit f l <C I2k, the final result is 


47 iRpu' 


(6) 


where o = dfl/dlnr is the shear rate in the WD. Note that 
this can also be rewritten in the dimensionless form 


a 

O 


fix 

n 


(7) 


This shows that a non-zero shear must exist throughout the 
star to allow the star to spin up and that it is larger when the 
viscosity is small (so that t v \ sc is big), the accretion rate is 
large (so that f acc is small), the thickness H of the layer is 
smaller, or when the spin 12 is smaller. 


Material accreting onto a WD at a rate M reaches its surface 
with a nearly Keplerian spin frequency 

n K = (GM/f? 3 ) 1 / 2 = 1.0 Ml /2 R~ 3 7 /2 s“\ (1) 

where Mi = M/1 Mq and i?g .7 = R/5 x 10 s cm are the typ¬ 
ical WD mass and radius I will be considering. Most of the 
energy of this incoming material is dissipated in a thin bound¬ 
ary layer and does not reach far into the star (Piro & Bildsten 
2004). Nevertheless, this material adds angular momentum 
at a rate of (GMf?) 1//2 M, and so a torque of this magnitude 
must be present to transport the angular momentum through 
the WD and increase its spin. 

Since these shallow surface layers of the WD have a typical 
pressure scale height H <C R, where H = P/pg and g = 
GM / R 2 is the surface gravitational acceleration, I assume a 
plane-parallel geometry and set R as the radius and use z as 
the vertical height coordinate above R. This also allows me 
to introduce a useful variable in the column depth 

dX = —pdz, (2) 


to measure the depth in the atmosphere. In the plane-parallel 
limit, the pressure given by hydrostatic equilibrium is sim¬ 
ply P = Xp. In this geometry, the transport of angular mo¬ 
mentum is described by a one-dimensional diffusion equation 
(Fujimoto 1993) 


d 

dt 


(R 2 n) 


1 d 
~R 2 p~dz 



(3) 


where 12 is the spin of the WD and v is the viscosity for trans¬ 
porting angular momentum. In principle v could be some¬ 
thing like a molecular diffusivity, but this is never relevant for 


2.2. Turbulent Mixing Framework 

Traditionally mixing is modeled as a diffusive process 
where for some element i with mass fraction X,, 


dXj = J_d_ f 2 dXA 
dt r 2 p dr \ ^ dr ) ’ 


( 8 ) 


where p is the density and D is the diffusion coefficient for 
mixing, which will vary with depth in the star. While this 
is well-suited for implementation in one-dimensional stellar 
evolution codes (e.g., Heger et al. 2000), an alternative way 
to think about the mixing is like convection. Just as a region 
in the star is either convectively mixed or not, a region in the 
star can be turbulently mixed or not. 

The general picture I employ is summarized in Figure 1 and 
described next. After accreting at a rate M for a timescale 
f acc , the column depth at the base of the accretion is 


X 


acc 


Mt a 


AnR 2 


(9) 


as shown by the horizontal dashed line in Figure 1. The tur¬ 
bulent mixing time at this depth is estimated as 


fmix(X acc ) = H 2 /D , (10) 

where the right-hand side is evaluated at X acc . If 
imix(X acc ) > f acc , then the mixing takes too long to act in 
comparison to the accretion and there is minimal mixing. If 
instead f m i x (X acc ) < £ acc , then the mixing can carry material 
to depths larger than the accretion depth. 

Evaluating f m ; x at every depth in the star, one can then find 
the depth where f m j x = f acc , which is denoted X m ; x . This is 
the depth down to where material is mixed, as shown in light 
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Figure 1. Diagram summarizing the main features of mixing considered here. Material is accreted for a timescale face at a rate M , which corresponds to a 
column depth E aC c = Mtacc/47r/? 2 . As long as f m ; x < face at E a cc. material is mixed down to the depth where f m ; x = face. This corresponds to a column 
of E m ; x > Eacc- The surface layer becomes completely mixed for all depths shallower than this (much like throughout the convective region of a star), and 
thus it has a helium mass fraction of T' m j x = YoE ac c/E m i x , where Yq is the helium mass fraction in the freshly accreted material. Carbon and oxygen are also 
mixed up into the surface layer and have corresponding mass fractions as specified in the diagram. 


gray in Figure 1. I assume everything is uniformly mixed 
down to this depth, which is roughly correct because the mix¬ 
ing timescale f m ; x is generally smaller at shallower depths. 
Thus if something mixes just a little at the base of the mixed 
region, the mixing is fast enough to mix it everywhere above 
this depth. In a sense, this is analogous to the widely used 
assumption of complete mixing in the convective regions of 
stars. Since material is mixed down to E ttlix , the mass frac¬ 
tion of helium in the mixed layer is decreased to be 


Unix — Yi E acc / Y 


(ID 


where Yq is the initial helium fraction of the material when it 
is first accreted (this variable is introduced just for complete¬ 
ness, I always use Yo = !)• Similarly, the C/O, which have 
mass fractions of Xy 2 and X\ (j in the C/O core, will be mixed 
up and have mass fractions 

X\2 ,mix — (1 — ^ acc /Emix), (12) 


like from Eddington-Sweet circulation). I hold off on con¬ 
sidering magnetohydrodynamic instabilities for future work. 
This would likely be something like the Tayler-Spruit dynamo 
(Spruit 2002), although this may overestimate the turbulent 
viscosity in some cases (e.g., Cantiello et al. 2014). 


3.1. Kelvin-Helmholtz Instability 

The Kelvin-Helmholtz instability (also referred to as the dy¬ 
namical shear instability) is governed by the Richardson num¬ 
ber 


Ri = 


N 2 


where N is the Brunt-Vaisala frequency 


N 2 


g_xr_ 

Hx P 


V a( J 


UhlT\ ' 

VdiZpJJ ’ 


(14) 


(15) 


and 


2fl6,mix — S acc /E m i x ). 


in the mixed layer. 


where \Q = 9 In P/d In Q, with all other intensive variables 
set constant, V a d = (9 In T/9 In P) a d is the adiabatic tem- 
(13) perature gradient, and the star refers to derivatives of the en¬ 
velope’s profile. This is just the thermal contribution to the 
buoyancy, and it is estimated to be (Bildsten 1998) 


3. TURBULENT TRANSPORT MECHANISMS 

To investigate the angular momentum and material trans¬ 
port discussed above requires setting v and D, respectively. 
For the present work I focus on transport mediated by hydro- 
dynamic instabilities triggered by accretion at the WD sur¬ 
face. I next summarize the Kelvin-Helmholtz instability and 
the baroclinic instability since these traditionally are most 
effective at radial transport (rather than latitudinal transport 


w ' s (ll) 1,2 = 2fiT2l/2s - 1 - (16) 

Note that for the scalings presented in this section I focus 
on the thermodynamic properties of the surface layers and 
omit the scalings with the WD mass, radius, and composi¬ 
tion for simplicity. These scalings all assume M = 1 M 0 , 
R = 5 x 10 s cm, and a helium-rich composition, but I 
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consider changes from these WD properties in the numer¬ 
ical investigation later. Linear analysis shows that Kelvin- 
Helmholtz instability occurs when Ri < 1/4, which develops 
into strong turbulence that readily transports angular momen¬ 
tum. 

The condition for instability of Ri < 1/4 assumes that 
thermal diffusion can be ignored for the unstable fluid per¬ 
turbations, in other words, that the perturbations are adia¬ 
batic. Fluid perturbations with a characteristic size L and 
speed V become non-adiabatic when the timescale for ther¬ 
mal diffusion, L 2 / K, where K is the thermal diffusivity, is 
less than the timescale of the perturbation, L/V. The ratio 
of these two timescales is the Pec let number, Pe = VL/K 
(Townsend 1958). The restoring force provided by thermal 
buoyancy is weakened when Pe < 1, which requires the sub¬ 
stitution of N 2 —>• PeN 2 and promotes instability. Thermal 
diffusion is most efficient at small lengthscales, which mo¬ 
tivates setting LV/v k = Re c (Zahn 1992), where v k is the 
kinematic viscosity and Re c is the critical Reynolds num¬ 
ber for turbulence, which is of order 1000. This gives the 
Peclet number approximately related to the Prandtl number, 
Pr, by Pe « Re c Pr. The turbulent perturbations are thus 
non-adiabatic when (Zahn 1992) 

K > v k Re c . (17) 


In the non-degenerate surface layers the kinematic viscosity 
is dominated by ions, and has a value of (Spitzer 1962) 

v k = 1.4 x 10 _2 pg ^g 572 cm 2 s" 1 , (18) 

where p$ = p/10 5 g cm -3 , and I assume a Coulomb loga¬ 
rithm of In A = 20. Setting K = IGctsbT ' 3 / (3c p np 2 ), where 
USB the Stefan-Boltzmann constant, c p the specific heat, and 
k the opacity, the thermal diffusivity is 

K = 980/C(/^ 2 T 8 3 cm 2 s"\ (19) 

where I approximate c p = T)k n j'lp:<n v and scale the opac¬ 
ity to rto .2 = k/ 0.2 cm 2 g _1 , appropriate for electron scat¬ 
tering in hydrogen-deficient material. Substituting Equations 
(18) and (19) into Equation (17), I find that the perturbations 
are non-adiabatic at depths of p < 7 x 10 6 g cm " 3 T s 1/2 . The 
new “secular” Richardson number associated with this limit 
is, 


Ri s 


v k Re c N 2 
K a 2 ’ 


( 20 ) 


When Ri s < 1/4, the so-called “secular shear instability” 
arises. 

The competing effects of accretion increasing a versus tur¬ 
bulence developing when Ri s <1/4 (and decreasing a) drive 
the surface layers toward marginally satisfying Ri s = 1/4 
(assuming for the moment that the sole viscous mechanism is 
the Kelvin-Helmholtz instability). This expectation is borne 
out in the white dwarf studies of Yoon & Langer (2004b). 
Thus I estimate the q due to this mechanism. Substituting 
Ris = 1/4 into Equation (20), and assuming Re c = 1000, 
the shear rate is 

_ on 1/2 l/2„—3/4 1 

tTKH — 2.0 a? 0 2 /?5 Tg S . (21) 


A shear rate this large is similar to the Keplerian frequency 
given by Equation (1), but as I show next, other fluid instabil¬ 
ities limit the shear before it can reach such a large value. 


3.2. Baroclinic Instability 

The baroclinic instability arises because surfaces of con¬ 
stant pressure and density no longer coincide if hydrostatic 
balance is to be maintained when differential rotation is 
present. In such a configuration, fluid perturbations along 
nearly horizontal directions are unstable, though with a suffi¬ 
cient radial component to allow mixing of angular momentum 
and material. The instability can roughly be broken into two 
limits, depending on a critical baroclinic Richardson number 
(Fujimoto 1987), 


Ri bc = 4 




lOJTafl 


o.l) 


( 22 ) 


where fio.i = £2/0.1 s -1 . When Ri > Rise, Coriolis effects 
limit the horizontal scale of perturbations. This results in two 
parameterizations for viscosity estimated from linear theory 
(Fujimoto 1993), 


A'BC 


«BC 

~ 3 ~ 


1 


H 2 n, 


a bc Ri bc 
3 Ri 3 / 2 


H 2 fl, 


Ri < RiBCi 
Ri > 


(23) 


where I include a dimensionless factor ubc, to account for 
uncertainty in how linear theory relates to the saturated ampli¬ 
tudes of the instability. In general, I find «bc 1 is required 
for significant mixing, which I explain in more detail below. 

By substituting vbc into Equation ( 6 ), I solve for the shear¬ 
ing profile. Due to the relatively low RiBC for WDs, see 
Equation (22), it is almost always the case that Ri > Pine- 
Thus, using the second of the two viscosity prescriptions I 
find 

Ri = 3.6 x 10 4 a^ /2 p^ /2 T 8 1/4 Q 3/ 1 2 MI 6 1/2 , (24) 


where M _ 6 = M/10 6 M 0 yr l . This demonstrates di¬ 
rectly that Ri 3> Ribc for accreting WDs. The shear is 


_ -i <„in-2 -1/4 5/8 n —3/4 ,Vl/4 _i 

atbc — 1-4 x 10 cr BC p 5 7 g flo.i ^f -6 s > 


(25) 


This demonstrates a couple of key features of the shear. 
First, it is higher for large accretion rates, which makes per¬ 
fect sense because angular momentum is being added more 
quickly to the star. Second, it is smaller for higher spin rates. 
This is because the baroclinic instability becomes stronger 
when the WD is spinning faster, which leads to a smoothing 
out of the shearing. Also note that generally ctbc ctkh, 
so that the baroclinic instability triggers before the Kelvin- 
Helmholtz instability. This prevents the shear rate from ever 
becoming large enough for the Kelvin-Helmholtz instability 
to operate at depths of p > 100 g cm -3 , i.e., at depths critical 
for unstable helium ignition. 


3.3. Viscous Heating 

The shear and dissipation by viscosity also leads to viscous 
heating within the accreted layer. The heating rate per unit 
mass is 

e = vo 2 / 2 , (26) 

which for the shearing estimates given above results in 

e B c = 6.9 x 10 2 a- 4 / 4 p- 5 / 4 T 8 - 5 / 8 £2-? / 4 M! / 6 4 ergg- 1 s- 1 . 

(27) 
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To put this in context, it is helpful to multiply this result by 
the accretion timescale t acc to roughly give the total energy 
per nucleon accreted 

E bc » 0.04a“c /4 P5 1/4 r 8 3/8 f2o.? /4 ^-6 4k eVnucr 1 . 

(28) 

In comparison, the thermal energy at the base of the accreted 
layer is about « lOkeVnucl” 1 . This means that viscous 
heating is usually not important, but can be for especially 
low values of a B c and (as can be see in discussions by 
Yoon & Langer 2004a, but in the context of shear instabili¬ 
ties). For the remainder of this work I ignore viscous heating 
so as to focus on the impact of mixing, but heating should be 
included in a fuller, more self-consistent calculation. 

3.4. Turbulent Mixing Estimates 

The turbulence that transports angular momentum also 
mixes material, albeit with much less efficiency because it re¬ 
quires more work to exchange fluid elements than to just exert 
stresses. Energy arguments give a mixing diffusivity that is 
related to the viscosity by (Piro & Bildsten 2007) 

D « v/Ri, (29) 

where there is a correction of order unity in this relation (often 
referred to as the flux Richardson number Rf), which I ignore 
for simplicity. I already have one free parameter in a B c t° 
adjust the strength of the turbulence (although this does not 
allow me to vary the angular momentum and material mixing 
independently). 

For this mixing diffusivity, I can now evaluate the mixing 
time as a function of depth using Equation (10), resulting in 


model, which extends down to the depth helium burning ig¬ 
nited unstably. This condition is given by the relation 

d( 3r> dc coo \ / ^ \ 

~dT ~ cLT ’ ( ’ 


where e 3a is the triple-a heating rate (Fushiki & Lamb 1987) 
and 


4 t 7 SB T 4 

£co ° 1- 


(34) 


is the radiative cooling rate. Both derivatives are taken at con¬ 
stant pressure. I solve for the temperature profile by solving 
the radiative diffusion equation. 


16cr SB T 3 dT 
3k d£’ 


(35) 


where I set the flux to a value which results in the correct ig¬ 
nition depth for the given WD mass and accretion rate. Since 
I am ignoring viscous heating, I can use the same tempera¬ 
ture profile whether or not turbulent mixing is included. The 
opacity is dominated by electron scattering at the relatively 
hot temperatures of the high accretion rates I consider. 

When mixing is included, it is incorporated as described in 
Section 2.2. Summarizing here, I let the WD accrete for a 
time face. I then solve for the mixing depth with an iterative 
process until the base of the layer correctly satisfies i m ; x = 
i acc . Once I know the depth of the mixing, I can compare 
the base conditions to Equation (33), making sure to use the 
correctly found Vj tlix in the stability criterion. If deza/dT < 
de coo \/dT, then the layer is stable and accretion can continue. 
Once dez a /dT = de coo \/dT, then the layer unstably ignites, 
and I can record the amount of mixing in that given model. 


fmix = 2.1 x 10 4 a^ 4 ^ /4 T 8 13/8 O^ /4 M_- 6 5/4 yr. (30) 
In comparison, the accretion time down to a similar depth is 
f acc = 4 ttR 2 H P /M = 1.9 x 10 3 p 5 T 8 Mr 6 1 yr. (31) 

This shows that / m ; x and £ acc are somewhat comparable for 
reasonable parameters, motivating that mixing is potentially 
important. Furthermore, since f m i x /f a cc -M -1 / 4 , the mix¬ 
ing gets stronger as the accretion rate increases. Conversely, 
the scalings with density argue that mixing is only effective in 
the layers sufficiently shallow that f m ; x < £ acc or when 

p < 5.9aB 3 T 8 " 5/32 0^ 3 M_ 6 gem" 3 . (32) 

This relatively small number for p means that both the viscos¬ 
ity parameter must be small a B c ~ 10“ 2 and the spin must be 
relatively small 12 ~ 1 0 _ 2 12 k for mixing to happen. This mo¬ 
tivates the values I use for these parameters for the numerical 
calculations in the next section. In particular, the large expo¬ 
nent of 12 -3 shows that the spin plays an especially important 
role in determining whether mixing is occurring. Also, note 
that the mixing I find here is never as large as what Yoon et al. 
(2004) calculate when focusing on the secular shear instabil¬ 
ity. This means that I also do not find solutions where the 
helium is mixed sufficiently enough to prevent unstable shell 
ignition. 

4. NUMERICAL CALCULATIONS 

I next explore the mixing in more detail with numerical 
models. For purposes of the present study, the profile for these 
conditions can be estimated with a constant flux atmosphere 


4.1. Example Accumulating Model 

To make this scheme more concrete, it is helpful to to con¬ 
sider a specific example of the time evolution of an envelope 
as it accretes toward unstable ignition. This example is for a 
M = 1 Mq WD accreting at a rate M = 10 -6 Mq yr -1 . A 
reasonable ignition mass for such conditions is Mj gn ps 3.0 x 
10~ 3 Mq (Iben & Tutukov 1989; Shen & Bildsten 2009). For 
the mixing, I set the WD spin to be 12 = 0.01 Ok and use 
cn B c = 0 . 01 . 

Figure 2 shows four different moments in time as the WD 
accretes and builds a surface layer of accreted material mixed 
with C/O from the WD surface. In each panel, the amount 
of time for which the WD has been accreting t acc is denoted 
along with the total mass fraction of C/O mixed into the sur¬ 
face layer. The solid line shows the temperature profile down 
to the base of the layer, with the solid circle denoting the 
amount of material that was actually accreted (much like the 
horizontal dashed line in Figure 1). The dashed curve is the 
stability criterion given by Equation (33). In the bottom right 
panel the surface temperature profile now extends down to this 
ignition curve demonstrating that this model will now ignite. 

A number of important features and trends caused by the 
mixing, that were described previously in my analytic explo¬ 
ration in Section 2, are now made more clear by this numerical 
example. These are as follows. 

1. Mixing is strongest at earlier times and decreases the 
longer the accretion persists until ignition (as can be 
seen by the mixed C/O mass fraction decreasing with 
time). This because it gets more and more difficult to 
mix to larger depths. 
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Figure 2. Four snapshots in time as an accreted layer grows toward unstable ignition. For this example, M = 1 Mq, R = 5 X 10 s cm, il = 0.01 Ok ■ and 
a BC = 0.01. Each panel shows the envelope temperature T versus column depth S profile as a solid line. The dashed curves indicate where unstable ignition 
occurs, which is found using Equation (33). The filled circles indicate the depth to which the envelope would extend if mixing was not included. The total mass 
fraction of C/O material mixed into the surface helium layer is denoted in each panel. Note that I denote the mass fraction variables as X 12 and A'i6 even though 
these refer to the mixed fraction (e.g., I have removed the subscript “mix” for brevity). 


2. Mixing causes unstable ignition of the surface layer in 
a larger mass of material for a given accretion rate than 
without mixing. In this case the mass of the surface 
layer is 3.4 x 10 -3 Mq at the moment of ignition after 
accreting 2.1 x 10~ 3 M Q of helium. This is « 13% 
more mass in the surface layer than the unmixed case. 

3. Igniting deeper also decreases the time to accrete until 
ignition, in this case by « 30% from « 3 x 10 3 yrs 
down to ~ 2.1 x 10 3 yrs. 

Most importantly of all, this demonstrates that it is reason¬ 
able for ignition to occur in a fairly mixed environment. In 
this case, Y ]t ,\ y = 0.62 while the mass fraction of C/O is 
0.38. This hopefully motivates the investigation of simi¬ 
larly mixed conditions in future surface burning models (e.g., 
Shen & Moore 2014). 

4.2. Trends with Accretion Rate and Spin 

To further illustrate the trends of the mixing, I explore the 
amount of mixing as a function of the accretion rate and WD 
spin. In Figure 3, I use M = 1 Mq, R = 5 x 10 s cm, and 
ctBC = 0.01 and plot contours of the mass fraction of C/O. 
This shows there is quite a range of mixing, from nearly un¬ 
mixed on the right side of the plot to fully 90% C/O on the 
left side. 

The general trend is that more mixing occurs at high accre¬ 
tion rates and low WD spin. The impact of a high accretion 
is intuitively clear, because if more angular momentum and 
material is added more quickly it is reasonable to expect more 



n/o K 

Figure 3. Solid lines are contours of constant mass fraction of C/O in the 
mixed accreted layer as a function of the accretion WD spin O and the ac¬ 
cretion rate M. This calculation uses M = 1 Mq, R = 5 x 10 8 cm, and 
ubc — 0.01. These contours are spaced every 0.1 in mass fraction as la¬ 
beled. The dashed curve shows the boundary at which to the right less than 
0.02 of C/O is mixed up into the accreted material. The wiggles in the dashed 
line are numerical and do not represent a real physical effect. 

mixing. Less clear is why slower spin leads to more mixing. 
The first thing to note is that it is not because there is a larger 
shear between the WD and the accreted material; the total 
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n/n K 


Figure 4. The same as Figure 3, but with a = 0.1. 

change in spin is basically the same whether fl = 0.00117k 
or 0.0117k- It is instead because the viscosity scales with the 
spin as shown in Equation (23). As the WD spins higher, sur¬ 
faces of constant pressure and density become more and more 
misaligned, which in turn drives the baroclinic instability and 
turbulent viscosity associated with it. A higher viscosity re¬ 
sults in a smaller shear to transport the angular momentum, 
as shown in Equation (6), and less shear leads to less mixing. 
Conversely, for small spin and turbulent viscosity, the shear is 
large and thus the mixing with it. This is particularly striking 
in Figure 3, where to the right of the dashed line there is quite 
little mixing. 

Since the analytic estimates showed a strong dependence 
on aec in Equation (32), I present the mixed C/O mass frac¬ 
tions for ckBC = 0.1 in Figure 4. This roughly moves all the 
curves to the left and up, showing how much more difficult it 
is to mix when the viscosity is higher. This is simply because 
there is less shear and turbulence in this case. For anc ~ 1 
there is little mixing except in the most extreme circumstances 
(of high accretion rate and low spin), and thus it is possible 
that turbulent mixing is negligible depending on what the real 
turbulent viscosity should be. 

As mentioned in the example calculation above, the mix¬ 
ing also causes ignition deeper than without mixing. This is 
quantified further in Figure 5, where I plot contours of con¬ 
stant ignition mass. This shows that at fixed accretion rate, 
the ignition mass goes up for lower spin (i.e., moving to the 
left). At fixed spin, the ignition mass goes down for higher 
accretion rate just as in the non-mixed case, but just not by as 
much. 

4.3. Trends with White Dwarf Mass 

The strength of the mixing should also depend on the mass 
and radius of the WD. To explore this, I repeat my analysis of 
for the amount of mixing which is summarized in Figure 6. 
This shows that there is in fact a strong dependence on mass. 
This can be understood from two main factors. First, the Kep- 
lerian spin of the incoming material is larger with larger mass, 
which means a larger amount of angular momentum needs to 
be carried through the star. Second, the larger surface gravity 
in turn results in a smaller scale height H in the surface lay- 



n/n K 


Figure 5. The same as Figure 3, but the contours are instead for constant 
ignition mass. Each curve is spaced by 2 X 10“ 3 Mg, some of which are 
labeled. 



Figure 6. Similar to Figure 3, but now mixed C/O mass fraction as a function 
of the WD mass. 


ers. This weakened the turbulent viscosity, leading to a larger 
shear and more mixing. 

5. DISCUSSION AND CONCFUSIONS 

I have explored mixing in the surface layers of helium- 
accreting WDs via hydrodynamic instabilities. These models 
highlight a number of important trends for mixing, which are 
as follows. 

1. Mixing is greatest at low spin f7, high accretion rate M, 
and high WD mass M. 

2. Mixing is stronger for low ctbc> an( I f° r ctBC ~ 1 mix¬ 
ing is negligible. 

3. Mixing causes unstable ignition of the surface layer to 
occur sooner and involve more mass. 
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4. Depending on the parameters, ignition of the surface 
layer can occur in highly mixed conditions with C/O 
mass fractions sometimes > 0.5. 

These results will hopefully motivate studies that follow the 
evolution of the surface layers after ignition occurs (e.g., 
Shen & Moore 2014), but including strong mixing and larger 
ignition depth. In particular, it will be important to under¬ 
stand how the nucleosynthetic products are different. Histor¬ 
ically, the double detonation scenario has had trouble with 
a high abundance of iron-peak elements near the surface of 
the star which causes the colors of the resulting SN la to 
not match those observed (Woosley & Kasen 2011; Sim et al. 
2012 ), but a large amount of mixing may alleviate this prob¬ 
lem. The different composition may also imprint itself on 
the properties of the rising light SN la light curve (Piro 
2012; Piro & Nakar 2013, 2014). Furthermore, mixing tends 
to make more calcium-rich burning products. In particular, 
Waldman et al. (2011) find that « 30% carbon mixed into the 
accreted helium may help explain the calcium-rich composi¬ 
tion in a subset of rapid transients (Kasliwal et al. 2012). 

The study presented here is mainly focused on highlighting 
the main properties of turbulent mixing, and thus there are 
many additional features that still need to be investigated in 
more detail. Most notably, the ignition conditions and enve¬ 
lope profiles were not calculated from first principles. Instead 
the known ignition conditions from the literature were used 
for the non-mixed cases, and then these were altered in a rea¬ 
sonable way to infer the changes from mixing. This can be 
improved by implementing mixing into accreting WD models 
in a stellar evolution code like MESA (Paxton et al. 2013), so 
that one can follow the evolution of the WD up to the point 
of ignition. The mixing would likely have to be implemented 
as a diffusive scheme (e.g., Heger et al. 2000). Although, as 
I argued earlier, from a physical perspective this is not nec¬ 
essarily a more accurate description of what occurs during 
turbulent mixing. Nevertheless, comparing this scheme to the 
one I use here is also an important exercise (such a compar¬ 
ison has been done in the case of accreting neutron stars by 
Piro & Bildsten 2007 and Keek et al. 2009, which shows at 
least qualitatively similar results). These studies could also 
explore the impact of other burning chains up until the point 
of ignition (besides the triple-a I focus on here) and viscous 
heating (see the discussion in Section 3.3). 

Another potentially important extension of this work is 
hydrogen-rich accretion and the subsequent classical novae 
(CNe) that result when mixing is included. Observation in¬ 
dicate that some ejecta of CNe are enriched by « 30% in 
C, N, O, and Ne (Livio & Truran 1994; Gehrz et al. 1998). 
Furthermore, a similar amount of enrichment has been ar¬ 
gued on theoretical grounds to explain the energetics of 
the ejecta (Starrfield et al. 1972; Truran 1982) as well give 
the best match with one-dimensional models of the ignition 
(Hernanzetal. 1996; Starrfield et al. 1998). The method for 
enriching CNe ejecta with heavier elements remains an out¬ 
standing theoretical question. If the enrichment is indeed im¬ 
portant for ignition, then it must take place prior to or during 
the early stages of the thermonuclear runaway. Previous the¬ 
oretical investigations include diffusion of hydrogen at low 
accretion rates (< 10~ 10 M 0 yr -1 , Prialnik & Kovetz 1984; 
Kovetz & Prialnik 1985). Others have looked at convective 
overshoot (Denissenkov et al. 2014), motivated by the work 
of Casanova et al. (2011). Turbulent mixing from accretion, 
as presented here, should also be explored as a possible solu¬ 


tion, and may be naturally stronger if accreting WDs are spin¬ 
ning slowly (Livio & Pringle 1998) or on ONe WDs (Mason 
201 1 ) because of their higher masses. 


I thank Ken Shen for useful exchanges on the double deto¬ 
nation scenario and surface burning on WDs, and I thank Ash¬ 
ley Ruiter and Stuart Sim for feedback on a previous draft. I 
also thank Mark Phillips for his generosity in arranging my 
visit to Las Campanas Observatory where much of this work 
was completed. 
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